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A MULTIPLE  RADAR  INTEGRATED  TRACKING  (MERIT)  PROGRAM 


INTRODUCTION 

During  fleet  exercises  in  the  1960’s,  it  was  demonstrated  that  many  targets  were 
never  detected  by  radar  operators  even  though  posttest  analysis  of  video  recordings  of  the 
radar  data  revealed  that  the  radar  return  from  targets  was  present  in  the  raw  video.  Among 
the  reasons  operators  missed  targets  were  operator  fatigue,  collapsing  of  upper  beams  of 
the  3D  radar  onto  a plan-position-indicator  (PPI)  display,  and  rain,  land,  and  sea  clutter. 

To  improve  its  surveillance  capability,  the  Navy  decided  to  associate  automatic  detection 
and  tracking  (ADT)  systems  with  its  radars.  Specifically,  the  SPS-48C  and  the  radar  video 
processor  (RVP)  for  2D  radars  have  been  approved  for  fleet  operation. 

On  board  most  naval  combat  vessels  are  two  kinds  of  surveillance  radars:  2D  and  3D. 
To  better  use  the  information  aboard  the  vessel,  the  radar  data  from  different  radars 
should  be  integrated  to  yield  a single  track  file.  The  benefits  of  such  a system  would  be 
increased  data  rate,  quicker  track  initiations,  increased  track  life,  and  on-line  redundance. 
Therefore,  in  the  spring  of  1973,  NRL  started  work  on  an  automatic  detection  and  inte- 
grated tracking  (ADIT)  system  [1-5]  that  performs  radar  integration.  Automatic  detectors 
are  associated  with  the  2D  and  3D  radars,  and  the  centroided  detections  are  integrated  into 
a single  track  file  in  a minicomputer. 

The  Applied  Physics  Laboratory  (APL)  embarked  on  a program  which  resulted  in  the 
SYS-l-D  system  [6].  The  SYS-l-D  system  adopted  the  NRL  system  philosophy  (merging 
detections)  and  uses  several  other  features  that  were  generated  by  the  ADIT  program.  The 
NRL  ADIT  system  is  a 6.2  research  program,  but  the  APL  SYS-l-D  system  is  a 6.4  program 
designed  for  fleet  introduction.  The  SYS-l-D  system  will  be  tested  in  the  fleet  during  1978, 

In  this  report  a large  computer  implementation  of  an  integrating  tracking  system  is 
discussed.  This  tracker  operates  on  suitably  formatted  detections  from  one  to  five  sensors 
yielding  a set  of  tracks.  The  tracker  is  essentially  the  ADIT  tracking  system  (4,  5]  . The 
motivation  for  this  implementation  is  twofold:  to  evaluate  radars  from  the  viewpoint  of 
their  performance  in  an  automatic-detection  automatic-tracking  mode,  and  to  improve  the 
performance  of  automatic  trackers  by  investigating  new  techniques  of  track  initiation, 
correlation,  and  filtering. 

Being  developed  concurrently  with  the  MERIT  system  is  the  SURDET  [7]  radar 
simulation.  This  program  simulates  the  detector  output  of  a radar  and  is  the  prime  source 
of  data  for  the  MERIT  program. 

An  overview  of  the  tracking  system  and  its  operation  is  given  along  with  an  outline 
of  possible  modifications  to  improve  tracking  under  various  conditions  of  input  data. 

BASIC  TRACKING  PHILOSOPHY 

Before  proceeding  to  the  principal  subject  of  this  report,  the  method  of  integration 
will  be  discussed  briefly.  Of  the  several  possible  methods  of  radar  integration  the  follow- 
ing ones  are  mentioned  most  frequently: 
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• Detections  from  both  radars  are  used  to  initiate  and  update  tracks; 

• Track  is  initiated  on  one  radar,  and  detections  from  both  radars  are  used  to 
update  tracks; 

• Tracks  from  both  radars  are  merged  to  yield  a single  track; 

• One  of  the  radars  is  chosen  using  a specified  criterion,  usually  the  highest  quality 
track  or  earliest  track. 


The  first  method  contains  the  maximum  information,  and  all  other  methods  can  be  con- 
sidered to  be  special  cases  of  it.  Consequently,  the  first  method  will  be  used. 

Three  types  of  tracks  are  considered;  clutter  points  (or  slowly  moving  targets),  target 
tracks,  and  tentative  (or  new)  tracks.  The  tracks  are  correlated  (associated)  with  the  detec- 
tions from  the  radars.  The  tracks  are  smoothed,  and  each  target’s  position  is  predicted  for 
the  next  time  one  of  the  radar  beams  will  pass  over  the  target.  To  reduce  the  number  of 
correlations  to  be  performed,  the  tracks  are  stored  in  azimuth  sectors,  and  only  those  detec- 
tions in  the  sector  where  the  track  is  located  and  in  neighboring  sectors  need  be  considered. 

Most  of  the  single-radar  tracking  systems  use  the  radar  itself  for  a clock,  since  the 
radar  operates  at  a constant  scanning  rate.  Although  this  system  is  similar  to  other  track- 
ing systems  using  a single  radar,  it  differs  from  previous  single-radar  tracking  systems  in 
timing,  filter  updating,  and  track  initiation  as  well  as  in  the  use  of  detections  from  several 
radars.  In  the  next  section  the  basic  system  parameters  are  defined  and  the  basic  routines 
are  discussed. 

TRACKING-SYSTEM  STORAGE  FILES  AND 
BASIC  ROUTINES 

When  a track  is  esUblished  in  the  software  of  the  computer,  it  is  convenient  to  assign 
a track  number  to  it.  With  this  system  all  parameters  associated  with  a given  track  are 
referenced  by  this  track  number.  Each  track  number  is  also  assigned  a sector  (region  of 
space  in  azimuth)  for  efficient  correlation.  In  addition  to  the  track  files  a clutter  map  is 
maintained.  A clutter  number  is  assigned  to  each  stationary  or  very  slowly  moving  target. 
All  parameters  associated  with  a clutter  point  are  referenced  by  this  clutter  number. 

Again,  each  clutter  number  is  assigned  a sector  in  azimuth  for  efficient  TOirelation. 

Tracking-Number  and  Clutter-Number  Files 

The  track-number  and  clutter-number  files  are  the  same  as  thosi'  described  by 
Richeson  of  APL  [S) . The  parameters  for  the  files  are  as  follows; 


NT 

NTL 

FULLT 

NEXTT 

LASTT 

LISTT  (256) 


Track  number. 

Last  track  number, 

Number  of  available  track  numbers. 

Next  track  number  available. 

Last  track  number  not  being  usetl. 

File  whosi’  256  locations  corre.spond  to  track  numbers. 
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NC  Clutter  number, 

NCL  Last  clutter  number, 

FULLC  Number  of  available  clutter  numbers, 

NEXTC  Next  clutter  number  available, 

LASTC  Last  clutter  number  not  being  used, 

LISTC  (256)  File  whose  256  locations  correspond  to  clutter  numbers. 

Only  the  operation  of  the  track-number  file  will  be  described,  since  the  operation  of  the 
clutter-number  file  is  identical. 


The  track-number  file  is  initiated  by  setting  LISTT  (I)  = 1 + 1 for  1 = 1 through  255. 
LISTT  (256)  is  set  equal  to  0 (denoting  the  last  available  track  number  in  the  file),  NEXTT 
is  set  equal  to  1 (the  next  available  track  number),  LASTT  is  set  equal  to  256  (the  last 
track  number  not  being  used),  and  FULLT  is  set  equal  to  255  (indication  that  255  track 
numbers  are  available). 

When  a new  track  number  is  requested,  the  system  checks  to  see  if  FULLT  is  0.  If 
FULLT  is  not  0,  the  new  track  is  assigned  the  next  available  track  number  (NT  = NEXTT). 
The  next  available  track  number  in  the  list  is  found  by  setting  NEXTT  equal  to  LISTT  (NT). 
FULLT  is  decremented,  indicating  that  one  less  track  number  ia  available.  Finally, 

LISTT  (NT)  is  set  equal  to  512  (a  number  larger  than  the  number  of  possible  tracks). 

This  is  not  necessary  but  helps  in  debugging  the  program. 

A track  number  is  dropped  by  setting  the  last  available  track  number  LISTT  (LASTT) 
equal  to  the  track  number  NT,  which  is  dropped.  LISTT  (NT)  is  set  equal  to  0 to  denote 
the  last  track  number,  and  LASTT  is  then  set  equal  to  the  track  number  being  dropj^ed 
(LASTT  = NT).  The  parameter  FULLT  is  incremented,  indicating  that  one  more  track 
number  is  available. 

The  track-number  and  clutter-number  files  maintain  a linkage  from  one  number  to  the 
next,  and  they  operate  rapidly,  eliminating  searching  techniques. 


Track  and  Clutter  Parameter  Files 


Parameters  associated  with  track  number  NT  are  as  follows; 


RTR (NT) 
TLT2CL (NT) 
RS  (NT) 

AS (NT) 

ES  (NT) 

VRS  (NT) 
VAS  (NT) 
VES  (NT) 


Smoothed  range  of  track,  stored  every  TTR2CL  seconds. 
Last  time  track  NT  checked  for  transfer  to  clutter  file. 
Smoothed  range  position 
Smoothed  azimuth  position  ^ JC*  (fe). 

Smoothed  height  position 
Smoothed  range  velocity 
Smoothed  azimuth  velocity  J Vg  (fe), 

Smoothed  height  velocity 


J 
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RPT  (NT) 

APT  (NT) 

EPT  (NT) 

TTR  (NT,  IRAD) 

TTL  (NT.  IRAD) 

TT  (NT) 

FA  (NT) 

FE  (NT) 

FVAR  (J,  NT) 


XMANT  (NT) 

TRKQAL  (NT) 
KT  (NT) 

TF  (NT) 


Predicted  range  position  | 

Predicted  azimuth  position  / 

Predicted  height  position  (•  ES  (NT)),  | 

Last  time  any  particular  radar  (radar  IRAD)  updated 
track  NT, 

Next  time  radar  IRAD  has  opportunity  to  uixlaU* 
track  NT, 

Last  time  any  radar  uixiated  track  NT, 

Azimuth  of  the  initial  detection  of  a tentative  track. 
Height  of  the  initial  detection  of  a tentative  track. 

Smoothed  average  of  sensor  measurement  errors  of 
sensor  uixlating  track  NT  in  range  (J  = 1),  azimuth 
(J  • 2),  and  height  (J  = 3), 

Time  since  the  initiation  of  tentative  track  NT  or 
the  time  since  the  last  maneuver  of  firm  track  NT , 

Time  since  the  last  3D  update  on  track  NT, 

Flag  indicating  a tentative  (-1)  or  firm  (0)  track 
status  for  track  NT, 

Time  of  initiation  of  track  NT. 


Parameters  associated  with  clutter  numlrer  NC  are  as  follows: 


RC  (NC) 

TLC2TR  (NC) 

RPC  (NC) 

APC  (NC) 

EPC  (NC) 

TC  (NC) 


Range  of  the  point  clutter  stored  every  TCL2TR 
seconds, 

Last  time  clutter  (roint  NC  was  checked  for  transfer 
to  the  track  file. 

Range  of  the  point  clutter. 

Azimuth  of  the  point  clutter. 

Height  of  the  point  clutter  (unused), 

Lasit  time  clutter  NC  was  updated. 


Tracking-Numlrer  Assignment  to  Azimuth-Sector  Files 


The  azimuth-range  plane  is  st'piuated  into  equal  azimuth  st'ctors.  After  a track  is 
utKlatixl  or  initiaU'd,  the  priniictwl  position  of  the  target  is  checked  to  see  which  sector  it 
occupies,  anil  the  track  is  assigiuHl  to  this  sector.  If  the  track  is  dropinxl  or  movi's  to  a 
new  si'ctor,  it  is  droppeil  out  of  the  sector  in  which  it  was  previously  locatixl.  The  para- 
mi'ters  associateil  with  sector  files  are  as  follows: 

ITBX  (I)  First  track  numU'r  in  sector  I (a  subscript  of  array  IDTl. 

IDT  (256)  Storage  location,  ixrrresponding  to  one  of  the  256  track 

numbers  and  containing  either  the  next  track  number  in 
.sector  I or  a 0. 


J. 
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U'BX  (1.  IKAl» 


IIH' 


First  rlutti'r  luinilH'r  in  wotor  I (u  sulisiTipt  of  urrtiy 
ini’)  from  nullur  IK. VO, 

Storago  loi-ation.  ixirrospomlmi;  to  oiu-  of  tlu-  'Jf'K 
i-hittor  mmilH'rs  aiui  lOiitmiimn  I'lthor  tiu'  iioxl 
I'luttiT  uuiiiIht  111  M'l-for  1 I'r  u 0. 


I'lio  assiKiimi'iil  Iraok  miiiilvrs  to  a/iiiiutli  ivlls  will  l>»>  il»>siril>»\l,  tiio  ilultor  iiiiiiitH'r 
assijti'ii'i’i't  IS  similar,  aiul  tlio  priuvss  is  os-si'iitially  tlu'  saiiu'  as  ilosi'nlnsl  iii  Kof.  »1.  ll'«' 

I TBXtl)  fill'  loiitaiiis  tlio  first  traik  luimlvr  iii  sivtor  I.  If  ri'HX(l)  - 0.  tlmn-  art*  no  tracks 
in  si-itor  I I'lu'  inr  (2r>t'>)  file  has  stora»;o  locations  cv>rrcspoiulin>;  {o  each  of  the  po.ssihle 
2f>t'  track  niniil>ers.  The  first  track  luimlx'r  m si'itor  I is  ohtaiiuHl  from  FIRST  “ lUlXtl). 
Till'  Ni'ixMul  track  niiinher  in  the  si'ctvir  is  I'htanuxl  hy  NFXTl  “ IPT  (URSl).  Ihe  next 
track  luimher  m the  sivtor  is  ohtain.xl  hy  NF.XT2  - UVTtNFXTn.  This  proivss  is  ixm 
tiiuieil  until  a 0 is  emx'iiiitereil,  iiulicatinn  that  there  are  no  nion>  track  iiiinilH'rs  in  the 
si'ctor. 

When  a new  track  is  aiKhxi  or  a track  moves  from  one  s«'ctor  to  another,  a track 
niimher  must  he  aiUleil  to  the  si'ctor.  When  the  track  is  a tentatiw  in«'''l  track,  the  track 
iiumlH-r  N’T  is  plaivil  on  the  hottoin  of  the  stack.  The  file  is  traixnl  as  ilescrilnxi  until  a 0 
IS  eiuxuinten-d.  If  this  occurs  m IIVT  (NKXT).  then  lOT  iNF.XT)  is  s*'t  equal  to  NT  aiul 
IIVT  tNT)  IS  s«'t  equal  to  0.  If  ther»'  were  no  tracks  m this  .si'clor,  then  ITBX  ll)  is  si't 
.xlual  to  NT  ami  ivaiii  IPT  (NT)  is  si't  .*qiial  to  0.  When  the  track  heitq;  tniiisferreil  into 
a si'ctor  IS  a firm  track,  the  first  track  niimher  in  the  si'ctor  is  ston'd,  the  track  numher  NT 
iH'iiiK  addixl  IS  made  the  first  track  niiiiilH’r  in  the  si'ctor,  and  tlie  track  niimln'r  NT  in  the 
IIVT  (N'T)  file  is  made  equal  to  the  original  first  track  nunilH'r  in  the  si'ctor.  This  dual 
(iroixnlun'  ensun's  that  first  the  firm  tracks  an*  (irocesst'd  fnnii  the  toi'  of  the  file  and  then 
the  tentative  tracks  an'  proix'ssi'd  fnmi  the  bottom  of  the  file. 

When  a track  is  ilro|>|>«'d  or  movixl  out  of  a s»'ctor  I,  the  track  numher  must  he  re- 
movixl  from  the  sector  files.  To  facilitate  this  n'lnoviil,  as  each  track  niimlH'r  N'T  is  pro 
ix'.s.si'il,  t)ie  last  track  nunilH'r  (irovx'ssml  is  savi'xl  in  N'Tl,.  If  NT  is  the  first  track  numlH'r 
in  the  m'ctor.  N'Tl,  is  0.  To  dro)'  track  numlH'r  N'T.  the  proivdun'  is  as  follows.  When 
N’Tl.  is  0.  TTBX  (I)  is  si'f  equal  to  IPT  (N'T),  and  IP'T  (N'T)  is  s»'t  equal  to  0.  When  N Tl, 
is  not  0.  IP'T  (N  Tl.)  is  si't  iHiuiil  to  IP'T  (N'T),  and  IPT  (N'T)  is  si'l  equal  to  0.  This  pm- 
vvss  shortcirciiits  the  linkiiKi's  throuKh  the  si'ctor  file,  eliminatum  the  track  numlH'r  N'T. 
St'ttmtt  IP'T  (NT)  equal  to  0 has  no  function  in  the  pn'pram  other  than  as  a dehunuinj!  aid. 
The  last  step  is  to  si't  N'T  equal  to  IP'T  (N'Tl.)  (or  to  TTBX  (P)  to  obtain  the  next  track  to 
Ih'  (inux'SjHHl. 

The  flowchart  for  m'ctor  file  mani|uilation  is  sliown  in  Fig.  1.  .X  track  does  not  change 
a'ctors  if  chaiiKinn  .-a'ctors  wouKI  ciiusi'  the  track  to  miss  an  op|Hirt unity  to  be  uiHlatiHl  hy 
another  nular. 

The  only  difren'iia'  U'twiH'n  the  ojH'ration  of  the  track-nomlH'r  assiKuments  and 
clutter-numiH'r  assignments  is  the  existamv  of  a clutter  file  for  each  railar.  'This  was  die- 
tatiHl  by  exiH'riments  on  the  .VPTT  system  Ifi)  usiiijs  the  .STS -I’i  and  the  STS  39  radars. 
Pbm'nation  of  the  detections  of  thesi'  two  radars  (one  on  l.-batul  radar  and  Ihe  other  on 
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S-batKl  radar)  demonstrated  that  few  of  the  f\xed  clutter  tH>int$  were  detected  hy  both 
radars.  The  first  clutter  numlK>r  in  each  sector  is  obtained  from  the  ICBX  file  ix^rres- 
punding  to  that  sector  and  the  radar  the  detection  comes  from. 


MSECT- CURRENT  SECTOR 
ISECT-SECTOR  IN  WHICH 
TRACK  BELONGS 
NT-CURRENT  TRACK 
NUMBER 

NTL-LAST  TRACK  NUMBER 
NE-NEXT  TRACK  NUMBER 
NL -TEMPORARY  VALUE 
NK- TEMPORARY  VALUE 


IS 

MSECT 

OUE  TO  BE  UPDATED 
BY  THE  OTHER 
RADAR 

(NO) 

I DROP 

rNFMOTlNT) 

I FROM  IDT  (NTL)  = NE 

I SECTOR 


NE=I0T(NT) 
ITBX(MSECT)--NE  I 

J 


NLOTBXdSECT) 
ITBX(1SECT)=  NT 
IOT(NT):NL 


_ ■ (NO) 

rNL  = ITBX(ISECT) 


NTL'NT 
NT  »IOT(NT) 
RETURN 


GET  NEXT  TRACK 
NUMBER 


NT«NE 

RETURN 


PLACE  TRACK  IN 
SECTOR  ISECT 


Fig.  1 - Soctor-filc  manipulation 


From  then  on,  the  operation  of  the  clutter  s»>ctor  files  is  identical  to  the  operation 
of  the  track  sector  files.  With  only  one  file  (IDC)  for  clutter  numlH*rs.  all  256  clutter 
numbers  are  available  for  use  no  matter  how  the  fixetl  clutU'r  iwints  are  distributeil  among 
the  radars. 


6 


NHL  KKPORT  H2t)0 


Input  OaUi  Bank 

riu'  first  st't  i)f  ilata  rt'aii  by  tlje  trarkiir):  systom  is  roail  only  onoo  aiul  tx>nsists  of  tiu' 
following  input  I'araiuetiTs: 

NRAl)  Number  of  raiiurs, 

RAIvrYP  (IKAD)  Typo  (21>  or  ai))  of  radar  IHAD, 

VHDT  (IRAO)  Rotatioir  rate  tradians/sl  of  radar  IRAP, 

IT  llRAH)  Initial  azimuth  sootor  of  radar  IRAH. 

Tbo  basic  input  data  from  the  radars  can  bo  brokon  into  two  oatonorios:  railar 
moasuromonts  and  control  paramot»‘rs.  Tbo  paranu'tors  assooialt'il  with  tho  input  ilata  art' 
as  follows: 


RM  (K.  IRAD) 

Range  measim'inent  of  the  Kth  detection  on  radar 
IRAn, 

AM  (K.  IRAH) 

.Azimuth  measurement  of  the  Kth  ileti'Ction  on 
nular  IRAI>. 

KM  (K.  IRADl 

Klevation  measun'ment  of  the  Kth  detection  on 
radar  IRAU, 

TM  (K.  IRAO) 

Time  of  the  Kth  detection  on  radar  IRAD, 

TMRK  (I,  IRAIH 

Time  radar  IR.AD  crv'sses  si'ctor  1, 

NP  (1.  IRAIH 

Position  of  the  (Kuiiter  in  the  ilata  buffer  of  radar 
IR.AD  in  si'ctur  1. 

NB  (1.  IRAIH 

Number  of  detections  from  radar  IR.AD  in  si'ctor  1, 

IT  (IRAIH 

('urrent  ,s«‘<'tor  of  radar  IR.AD. 

It  would  not  bo  noirsaary  to  retain  this  type  of  buffertnl  input  in  this  version  of  tho 
tracker,  liocaust'  tho  system  obtains  its  data  from  mass  stonifte.  However  it  is  amvenient 
to  retain  this  structure,  becaust'  the  logic  is  well  develofH'd  and  old  detections  are  available 
for  logic  analysis. 

The  program  expects  the  data  to  bo  made  available  from  one  .source,  with  the  data  , 

formatttxl  into  blocks  such  that  each  block  corresi>onils  to  a sector  of  data  from  a radar. 

Kach  block  bi'gins  with  a sot  of  ixmtrol  parameters,  the  radar,  the  sr'ctor,  tho  time,  the  ’ 

number  of  detections,  and  the  buffer  pointer.  Following  tho  ixmtrol  data  are  the  moasurt'-  t 

nionts  for  each  detection,  the  range,  the  azimuth,  the  elevation,  and  the  time  of  detection.  | 

t 

I'o  ensure  that  the  data  prt'st'nfiHl  to  MKRIT  are  in  the  proper  format,  a program  for  j 

merging  radar  data  files  (MRliRORl  n'ads  data  files  from  different  radars  (currently  radars  i 

simulattnl  by  SUROKT),  |>artitions  the  detections  into  st'ctors,  onters  (he  sectors  in  time.  < 

luul  generates  a data  file  with  the  pro(H'r  format  for  MKRIT.  This  program  ixnild  be  jl 

nuxlifitHl  to  operate  on  real  radar  data;  however  the  data  must  Im'  from  colocatixl  radars 

operatwl  concurn'iitly  to  be  meaningfvil.  i 
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TRACKING  SYSTEM 

The  tracking  system  has  bt*en  discussed  in  previous  reports  on  ADIT  (4,  5).  The 
m^or  change  in  the  system  is  to  . ' iw  data  from  one  to  five  radars  with  arbitrary  para- 
meters (pulselength,  beamwidth,  scan  rate,  etc.)  to  Ik*  inU‘gratt\i  rather  than  data  from  ‘ 
exactly  two  particular  radars.  Briefly,  the  nuyor  features  of  the  system  are:  a methml  oT 
aihusting  the  tracking  parameters  to  the  characteristics  of  each  radar,  an  individual  chitler 
map  for  each  radar  which  eliminates  deU'ctions  that  correlate  with  a fixetl  clutter  point, 
a tracking  filter  to  estimate  track  parameters,  a correlation  logic  to  associate  detections 
with  the  projw  track,  a method  of  initiating  new  tracks  in  the  system,  and  a method  of 
conveying  information  to  the  user  on  the  operation  of  the  system.  Each  of  the.s*'  topics 
is  covered  in  a following  subsection. 

Tracking  Program  Control  Parameters 

To  retain  as  much  flexibility  as  tnissible,  many  parameters  of  the  tracking  system  are 
defined  when  the  program  is  run.  The  parameters  are  generally  those  that  should  be  bastxl 
on  the  quality  of  the  detection  data.  The  parameters  and  their  explanations  are  as  follows. 

Miscellaneous  Parameters 

SCNFRT  (usually  Is).  If  a clutter  point  or  track  has  shiftixl  sectors,  the  same  radar 
could  attempt  to  u(xlate  it  twice  on  the  same  scan.  SCNFRT  is  a time  delay  us»M  to  in- 
hibit correlations,  and  would  change  only  in  exct'ptional  circumstances,  such  as  when  a 
radar  had  a very  fast  scan  rate. 

NSECT  (usually  64',.  This  is  the  number  of  azimuth  sectors  and  would  change  only 
for  a radar  with  such  pwr  azimuth  accuracy  that  a tletection  could  occur  nrore  than  one 
sector  away  than  the  sector  of  the  track. 

NSECSH  (usually  8).  The  rotation  rate  of  the  radar  is  calculated  and  smoothed 
every  NSECSH  sectors,  and  this  numlrer  of  sectors  would  change  only  if  some  feature  of 
the  radar  necessitateil  a change  in  frixjuency  of  update. 

ICBIAS  (usually  6).  This  is  the  lag  of  processing  behind  the  position  of  radar  and 
would  change  only  if  the  program  was  chai\ged  to  perform  additional  functions  such  as 
merging  deU'ctions  or  removing  biast's  between  radars. 

NKON  (usually  5).  This  is  the  total  number  of  track  wrrelation-rt'gior,  sizes,  which 
can  be  increased  to  10  to  allow  finer  discrimination  in  dense'  target  environments. 

IFLAG  (usually  0).  A value  of  1 for  IFLAG  forces  use'  of  all  NKON  correlation 
rc'gions.  A value  of  0 mc'ans  that  on  tracks  with  “good”  velocity  estimates  only  the 
smallt'r  correlation-region  sizes  art'  used. 

IPCTR.  This  is  the  output  control;  0 implies  firm  tracks  only,  1 implies  firm  tracks 
plus  firm-track  comments,  2 implies  firm  and  tc'ntative  tracks  only,  3 implies  firm  and 
U'ntative  tracks  plus  track  cximments,  and  4 implies  firm  and  tentative  tracks  plus  all 
cximments.  Track  information  consists  of  the  value's  of  all  the  critical  track  parameters 
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at  oach  opportunity  to  uixlate  tlu  vraik,  oomnuM\ts  aro  ronuirks  on  siKmfu-ant  ownts 
involvinn  a traik  or  oluttor  suili  as  In'inn  dropped,  bo\nn  i-onvertint  fix>in  a track  ti>  a 
(’liitt«'r,  or  )iotin>;  whu.•l^  detections  correlate  with  a clutter  point. 

Filler  Parameters 

FLTNlhM  (J.  IKAD)  and  Kl.TDNM  (J.  IKAI'O.  Thes.'  are  the  numerator  and 
denominator  list'd  in  on  eijuation  to  si't  the  filter  bandwidth  m the  ran(;e  (J  “ 1),  a/imutli^ 

(.1  2).  and  height  (J  * 3|  irackint;  filters  for  ra»iar  IKAD  as  a function  of  track  ••quality". 

\I  (usually  O.tU,  This  is  the  dampmn  factor  of  the  filter. 

DTKL  (usually  30  s).  This  is  the  maximum  time  between  elevation  updates  on  a 
track  for  j;ood  elevation  information.  The  pn'si'iit  value  may  be  alteriil  after  evaluation. 

lU'Hir  (usually  0.5  n.mi.l.  When  the  error  between  I'revhcted  ('osition  and  measured 
(Hisition  excet'ils  RTRIT,  a target  maneuver  is  indicati'il.  I'he  value  depeiuls  on  trackmi; 
accuracy. 

Track  "Qualtlv  ” Parameters 

-•  XMXNtNT  (usually  00  si.  This  is  the  largest  value  that  the  track  parameter  Xkt.VN  T 
is  allowiHl  to  achieve.  There  is  no  reason  for  this  inirameter  ti>  chanjte  with  currer.  ."itic. 

FIMANT  (usually  5.6  si.  A time  XMAN  T larser  than  the  time  FIM.\N  T indicates  a 
gooil  "quality  ’ tentative  track  and  fonvs  a nxluction  in  the  numlx'r  of  iurrelatiiu,  ri'nions 
used.  Usually  5.6  s represt'iits  the  minimum  time  since  initiation  to  obtain  (jood  velocity 
information. 

KMANT  (usually  10  s).  .A  time  XM.ANT  larger  than  the  time  KM.ANT  indicates  a 
hoik!  “quality”  firm  track  and  forces  reihiction  in  the  numln'r  of  correlation  regions  useil. 

TFMNT  (usually  2 s).  If  the  time  XMANT  is  less  than  the  time  'I'FMNT,  the  time 
liast'  is  too  short  for  useful  velocities  to  be  calciilatiHl  for  a tentative  track. 

FFMNT  (usually  0.8  s).  If  time  since  last  update  for  a firm  trick  is  less  than  FFMNT, 
the  timebase  is  too  short  for  a valid  correction  to  Ih'  made  to  velocity. 

Time  Parameters 

TCL'ZrR  (usually  20  si.  Kvery  TCL2TR  si'ininds  a clutter  is  check»\l  for  motion  by 
calculating  a range  w'locity  Ivased  on  the  current  range  and  the  stonnl  range  fivm  T('l,’21'R 
si'conds  ago.  Clutter  points  with  velocities  greater  than  lUFt'T  (IR.-Vni  will  In*  ixmvertixl 
to  tracks. 

TTR2CL  (usually  20  s).  This  is  the  time  l»as«>  for  calculating  a velocity  for  a firm 
track.  Tracks  with  velocities  smaller  than  DIFTC  will  1h'  wnwrliHl  to  clutter  points. 


TCMAX  (usually  30  s).  This  is  the  time  an  unupdated  clutter  point  is  retaimxl  and 
IS  a function  of  the  sean  rates  of  the  radars. 
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TTMAX  (usually  20  s).  This  is  the  time  an  unupdated  firm  track  is  retained  and  is 
a function  of  number  of  radars,  scan  rates,  probability  of  detection,  etc. 

TNMAX  (usually  15  s).  This  is  the  time  an  unupdated  tentative  track  is  reteined 
and  is  a function  of  number  of  radars,  scan  rates,  probability  of  detection,  etc.  TNMAX 
should  be  less  than  TFIX. 

TFIX  (usually  20  s).  An  update  to  a tentative  track  at  a time  subsequent  to  TFIX 
seconds  after  initial  detection  causes  an  assignment  of  the  track  to  the  clutter  map  or 
firm-track  status  based  on  a comparison  to  the  velocity  parameters  VRMIN  and  VAMIN. 
TFIX  should  exceed  TNMAX. 

TUPTD  (usually  12  s).  If  the  time  since  the  last  update  of  a firm  track  exceeds 
TUPTD,  then  more  correlation  regions  are  used  to  compensate  for  a less  reliable  predicted 
position.  TUPTD  depends  on  expected  error  in  predicted  position  and  the  effect  of  a 
target  maneuver. 

TICLUT.  From  the  time  the  tracking  system  starts  until  time  TICLUT , only  clutter 
points  are  established  to  avoid  an  initial  overload  of  the  tracking  system  with  tentative 
tracks  based  on  fixed  clutters.  The  elapsed  time  usually  equals  two  scans  of  the  slowest 
scanning  radar. 

Velocity  Parameters 

VRMIN  (usually  50  knots)  and  VAMIN  (usually  0.0001  rad/s).  Tentative  tracks 
with  range  velocity  less  than  VRMIN  and  azimuth  velocity  less  than  VAMIN  are  converted 
to  clutter  points;  otherwise  they  are  converted  to  firm  tracks, 

DIFTC  (usually  50  knots).  A track  whose  average  range  velocity  over  the  time 
interval  TTR2CL  is  less  than  DIFTC  may  be  converted  to  a clutter  point. 

DIFCT  (IRAD)  (usually  50  knots).  A clutter  point  maintained  by  radar  IRAD 
whose  average  velocity  over  the  time  interval  TCL2TR  is  greater  than  DIFCT  (IRAD)  is 
converted  to  a firm  track. 

Radar  Measurement  Variance  Parameters 

RVAR  (IRAD).  This  is  the  variance  in  range  measurements  on  radar  IRAD.  RVAR 
is  a function  of  compressed  pulse  length  P and  usually  equals  P^/12  expressed  in  nautical 
miles. 

AVAR  (IRAD).  This  is  the  variance  in  azimuth  measurements  on  radar  IRAD. 
AVAR  is  a function  of  antenna  beamwidth  B and  usually  equals  (B/10)2,  the  accuracy 
possible  with  the  signal-to-noise  ratio  necessary  for  detection. 

EVAR  (IRAD).  This  is  the  variance  in  elevation  measurements  on  radar  IRAD. 
EVAR  is  a function  of  antenna  elevation  beamwidth  E and  the  type  of  elevation  scanning 
employed.  As  a worst  case  (the  assumption  of  discrete  beam  po.sitions  with  a detection 
in  only  one  beam)  the  variance  equals  £2/12. 
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Clutter  Correlatuiii-Region  Size 

CRC  ([RAD).  This  is  the  range  wm^lotion-region  size  on  railar  IRAD.  Tlie  usual 
value  is  one  range  (.•ell  in  nautieal  miles. 

CAC  (IRAD).  This  is  the  azimuth  correlation-region  size  on  radar  IRAD.  There  is 
some  indication  that  this  is  more  a function  of  the  spatial  distribution  of  the  scatterers 
causing  the  clutter  detection  than  of  the  azimuth  lu-amwidth  of  the  radar. 

CEC  (IRAD).  This  is  tlie  maximum  elevation  a detection  fnim  radar  IRAD  may 
have  to  be  allowed  to  ujHlate  a clutU*r  point.  It  is  a function  of  elevation  accuracy  ami 
U'rrain  geometry. 

Track  Correlation  Sizes 

CRT  (1,  IRAD).  This  is  the  Ith  range  correlation-region  size  (n.mi.)  on  railar  IRAD. 

CAT  (I,  IRAD).  This  is  the  Ith  azimuth  correlation -region  size  (radians)  on  radar 
IRAD. 

CET  (1,  IRAD).  This  is  the  1th  elevation  correlation-region  size  (radians)  on  radar 
IRAD.  The  values  for  all  of  the  correlation-region  sizes  depend  on  the  radars,  the  environ- 
ment, the  scenarios,  and  the  operational  riKpiirements. 

Individual  Clutter  Maps 

An  individual  clutter  map  removes  from  the  iletection  files  of  the  associateil  railar 
any  deU'ctions  which  correlate  with  clutter  points  or  slowly  moving  targets.  The  decision 
to  implement  individual  clutter  ma|)s  was  Inised  on  ex|H'riments  with  the  ADIT  system 
iiwolving  an  L-band  and  an  S-banil  radar  jS).  Examination  of  the  data  files  revealed  that 
of  the  clutter  points  consistantly  detecteil,  only  a small  i)er<x>nU>ge  were  delect ihI  by  both 
radars.  Thus  a detection  from  radar  A will  not  update  a clutter  poitd  generatiHl  by  radar 
B and  will  be  available  to  uixlate  the  appropriate  track. 

The  mechanics  of  a single  individual  clutter  map  are  as  follows:  (liven  that  radar 
IRAD  is  due  to  update  .wctor  ID  of  the  track  map;  then  sector  ISECT  = ID  + 3 of  the 
clutter  map  of  radar  IRAD  .should  be  iiiHlatinl.  The  correlation  and  updating  is  performixl 
through  nesteil  iterations.  Eirst  each  eligible  clutter  point  is  obtainixt  and  a iurrelation 
is  attem|)ted  with  each  eligible  detection.  The  eligible  clutter  points  are  any  clutter  points 
in  sector  ISEtT  which  have  not  been  iRnlatixl  already  on  the  current  scan  of  radar  IR.AD, 
(This  ixruld  occur  when  a clutter  point  in  .sector  ISKCT  - 1 is  uiHlatiHl  with  a detection 
which  moves  it  to  si'ctor  ISECT.)  The  eligible  clutter  points  are  obtainwi  through  linkage 
provided  by  the  clutter  azimuth-si  ctor  files  (described  earlier).  The  eligible  detections  are 
obtainixt  from  the  input  data  files  in  the  si'ctors  1SK(T  - I,  ISEt’'!’,  and  ISKC*!  + 1.  Ihe 
detection  parameters  im*  locatixl  .siHiuentially  and  referemvd  u.sing  NT  (1,  RAD),  which 
is  the  last  detection  in  sector  I from  radar  IRAD,  and  NH  (I,  IRAD),  which  is  the  mmiber 
of  detections  in  .sector  I by  radar  IRAD.  Each  clutter  point  is  obtaineil,  checked  to 


11 


J.l).  WILSON 


dt'U'rinino  that  the  last  mnlate  was  at  least  one  sean  a»;o,  anil  then  retaineil  while  eaeli 
eli)tihle  detection  is  eoinfmrwl  in  ranne,  azimuth,  and  elevation  for  correlation.  Kacli 
detection  that  correlates  Ls  eliminattH.!  from  the  detection  file.  The  detection  nearest  the 
clutter  point  in  the  miLximum-likelihood  simse  is  ust'd  to  u|Hlate  the  clutter  usinn  a simple 
<v  filter  (n  “ 0.5,  the  aver»4!e  of  the  two  positions).  The  meaiun^  of  the  maximum-likeli- 
hood si'nse  is  that  the  error  metric  is  the  ranjte  error  sipiared,  dividiHl  by  the  variance  of 
range  measurement  error,  plus  the  azimuth  error  siiuiutihI,  divided  by  the  variance  of  the 
azimuth  measurement  error.  The  time  of  the  u|Hlate  is  savtxl,  and  the  si'ctor  of  the  clutter 
point  is  changtxl  if  necessary.  lhiu|H)atiHl  clutter  points  are  chcckinl  fi>r  elimination; 
uixlated  clutter  (Kiints  are  checked  for  motion  and  ixinversion  to  a track.  If  a clutter 
|K)int  from  the  clutter  map  from  radar  IR.Al)  is  converted  to  a track,  the  remaining  clutter 
maps  are  seiu'cheil  for  a ix)rresponiling  clutter  point,  which  is  dropped  to  keep  it  from 
stealing  detections  from  the  newly  creatiHl  track. 

In  summary,  the  clutter  files  store  the  locution  of  the  clutter  points  or  .slowly  moving 
targets  and  remove  detections  fn>m  the  ruilar  file  that  correlate  with  them.  If  too  long  a 
time  elapses  between  uiHlutes,  the  clutter  point  is  droppcil.  If  the  clutter  point  has  suffi- 
cient velocity,  it  is  cimverti'il  to  a track. 

t 

f 

Tracking  Filter 

The  tracking  filter  implementiHl  is  an  filter  with  a decreasing  bandwidth  and  a 
turn  detector.  The  iv-(l  filter  equations  usihI  for  range  and  azimuth  after  a firm  target  has 
lH*en  established  are 

* -V  1^*)  + ‘'X-  l.v„,  (X’)  - X,, 

c,  (fc)  - c,  (k  - 1)  +—  |.r„,  (X)  - .V,,  (A-)I. 

‘k 

i 

I 

.r,,  (k  -r  1)  - (fc)  + c,  tfc)  , 

where 


jf,  (fc)  ” smoothed  position, 

Xp  {k)  “ pivdicUxI  position, 

and  fl*  - system  gains, 

(Ar)  “ smoothed  velocity, 

Tf,  and  Ti,^^  _ sampling  periods,  and 

jf„,  (k)  “ measimnl  position. 


The  system  gains  and  li),  are  si't  by 

«fe 

ami 

0k 


2 - ttfc  ’ 
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when*  ^ is  tho  filter  damping  coefficient  and  to^  is  the  filter  bandwidth.  The  filter  band- 
width is  calculatwl  from 


_ FLTNUM  (J.  IRAQ)  FVAR  (J.  NT) 

“ FLTDNM  (J,  IRAD)  + XMANT  * VAR  (IRADl  ’ 

where  J » 1 or  2 refers  to  range  or  azimuth,  IRAI)  refers  to  updating  with  radar  IRAD, 
ami  NT  is  track  NT.  FLTNUM  and  FLTDNM  are  tracking  program  control  parameters 
associattxl  with  the  filter.  VAR  is  the  tracking  program  wntrol  parameter  for  the  mea- 
surement variant  in  range  or  azimuth  for  radar  IRAD.  FVAR  is  a past  history  of  the 
measurement  errors  of  the  radars  which  have  updated  track  NT.  FVAR  is  updated  by  a 
simple  a filter  with  the  measurement  variance  of  each  radar  as  it  updates  the  track. 

XMANT  is  the  time  interval  over  which  an  accurate  track  has  Iwen  maintained.  When  a 
track  is  made  firm,  XMANT  retains  the  time  since  the  track  was  initiated.  When  a track 
is  updated  with  a small  range  error,  XMANT  is  set  by 

XMANT  * MIN  (XMANT  + T^,  XMXMNT), 

and  when  a track  is  uinluted  with  a large  range  error  (indicating  a maneuver),  XMANT  is 
s*'t  to  0.11  s (an  arbitrary  small  time  for  debugging  purposes).  Specifically,  a range  error 
is  small  or  large  de()ending  on  whether  the  error  between  the  prerlictetl  range  and  the 
measurwl  range  is  less  than  the  tracking  prtrgram  control  parameter  RCRIT.  The  para- 
meU'rs  FLTNUM.  FLTDNM,  and  XMXMNT  are  chosim  to  yield  a range  of  bandwidth 
values  from  1/36  to  1/6.  Kx|H>riments  with  ADIT  indicate  that  coj,  may  be  too  small  for 
slowly  turning  targets;  that  is,  the  error  builds  slowly  until  the  bandwidth  is  suddenly 
incxt'aswl.  This  would  be  a matter  for  investigation  if  this  type  of  filter  is  riMaintnl. 

We  have  modifier!  the  filter  descr(bed  alwjve  to  operate  on  elevation  information. 

This  is  because  we  may  have  uixlates  on  both  2D  and  3D  railars  and  the  accurar'y  of 
elevation  measurements  is  poor  relative  to  the  range  of  elevations  we  expect  to  obsa'cve. 

If  the  ujKlate  information  is  derivrxl  from  a 2D  radar,  no  elevation  updating  is  possible. 

If  elevation  information  is  available,  the  smoothing  is  as  follows;  The  target  is  assumed 
to  be  constant  in  altitude,  and  the  filtering  is  done  in  height  not  elevation.  The  filter 
is  a variable-bandwidth  « filter.  1'he  velocity  is  as-sumixl  to  be  zero,  and  we  merely 
attempt  to  follow  the  track  height  using  the  « filter.  The  banrlwidth  is  a function  of 
range  rate  ami  elevation  angle  motlifieil  by  the  past  history  of  measurement-error  varinn<x*s 
dividetl  by  the  current  measurement-error  variance.  The  smoothing  gain  n is  calcvdattnl 

0*1-0  ” , 

where  TQ  is  the  time  sinw  the  last  uixlate  with  a 3D  radar.  With  large  TQ,  n 1 and 
the  smoothtHl  height  becomes  the  measured  height. 

Tentativt'  tracks  art'  uixlated  differently.  With  8t>  little  past  history  available  the 
principal  effort  is  to  obUiin  the  best  estimates  possible  for  the  track  for  u.w  at  such  time 
as  the  track  is  ma<le  firn).  To  acwmplisli  this,  the  initial  range  and  azinmth  arc  stonxl, 
at  each  U|xlate  the  smoolluxl  range  and  azimuth  are  si't  to  the  measured  range  and  azimuth. 
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and  tln'  moasured  heiKht  is  averaged  with  tin*  smoothed  lieiglit.  If  a sufficient  time  base 
exists  (time  since  initiation  ^treater  thiin  TFMNT),  the  range  rate  and  azimuth  rate  lU’e 
estimatrtt  using  the  storeii  values  at  initiation  anti  the  current  measureti  values. 

Track  Oirrelation  and  Updating 

Before  the  correlation  begins,  the  tracks  which  have  already  been  u|xlated  on  the 
current  scan  by  this  radar  are  flagged.  I’his  can  occur  when  a track  changes  st'ctors.  The 
tracks  are  obtained  through  the  operation  of  the  track  azimuth-sector  files. 

The  Uack  correlator  attempts  to  minimize  incorrect  correlations  by  multiple  pass»\s 
through  the  track  files  with  incrementwl  correlation-region  sizes.  The  track  is  updattxl  in 
only  one  s»'ctor  (primary  sector).  The  tracks  located  in  the  two  si’ctors  in  advance  of  the 
primary  st'ctor  (s»'condary  sectors)  are  correlated  with  detections  also  but  are  not  updatixl. 

If  a detection  correlates  with  a track  in  the  primary  sector  and  is  the  minimum-distance 
■orrelation  with  the  target,  the  detection  is  usixl  to  update  the  tmget  track  and  is  droppeii 
from  the  iletection  file.  A tentative  track  which  is  being  uixiatwl  is  transformeil  to  a tiir- 
get  track  TFIX  st'cond.s  after  initiation.  The  flag  denoting  that  it  was  a tentative  track 
is  removed.  If  the  track  is  in  a sc'condary  sector,  the  detection  is  flaggcxi.  This  flag  pro- 
hibits the  detection  from  updating  a track  in  the  primary  sector  on  some  later  pass 
through  the  track  file  using  a larger  correlation-region  size  on  the  basis  that  the  detection 
has  already  Ix'en  demonstrated  to  be  closer  to  some  other  track.  Each  track  which  cxirre- 
lates  with  a detection  is  also  flagged.  This  prevents  the  correlator  from  accessing  the  track 
again. 

The  finer  the  correlation-region  incremetits,  the  bedter  the  correlation  will  ix'rform 
in  terms  of  the  only  judgment  criterion  available  to  it;  the  proximity  of  detection  to  a 
prcxlicted  track  position.  In  most  instances  this  correlation  procxHlure  prevents  a track 
which  is  not  detectcxl  on  the  current  scan  from  stealing  a eletection  from  a nearby  tnu-k. 

When  more  than  one  detection  correlates  with  a track,  the  “clost'st”  tlelection  in  the 
maximum-likeliliood  senst*  is  useil  to  update  the  track.  The  definition  for  ma.ximum-like- 
lihood  is  similar  to  the  definition  given  in  the  discussion  of  the  clutter  maps  but,  if  pos.si- 
ble,  includes  a term  for  elevation  error.  This  is  possible  only  if  the  current  radar  is  a 311 
radar  and  if  the  track  has  been  updated  recently  (last  IITEL  seconds)  by  a 311  radar. 
The  coonlinates  of  this  detection  are  used  in  the  smoothing  filter  to  updat»>  the  track 
pjirameters.  Then  the  time  this  track  will  be  seen  by  the  next  radar  is  calculatetl,  and  the 
position  of  the  track  at  that  time  is  predicUxl. 

Before  a return  to  the  executive  the  tracks  must  be  prepartHl  for  the  next  correlation 
process.  First  all  flags  are  removeil  from  the  tracks,  anil  any  unflaggrxl  tracks  in  the  pri 
mary  sector  are  coasted,  using  the  stored  velocities,  or  droppeil  if  tot>  much  time  has 
elapstxl  since  the  lust  ujxlate.  Then  each  truck  that  has  I'hangixl  s»'ctors  is  placed  in  the 
ctrrrect  sr'ctor  file.  Firm  tracks  are  placiHl  at  the  top  of  the  sector  file,  tmd  tentative 
tracks  are  placet!  at  the  IxJttom.  This  insures  that  as  the  tracks  are  sixiuentially  processtxl 
by  the  correlator,  the  established  tracks  will  have  the  first  opportunity  t»)  be  uptlattxl  with 
each  incremental  correlation  region.  As  additional  effort  to  eliminate  falst'  corndations, 
especially  in  the  cast'  of  a misst'd  detection  on  a track,  restrictions  are  st't  tm  the  maximum 
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errors  allowed,  dependent  on  the  past 
number  of  correlation  regions  used  by 


history  of  the  track.  This  is  done  by  limiting  the 

each  track  through  the  following  logic:  i 


• A track  is  limited  to  using  four  regions  if  it  (a)  is  a tentative  track  with 
XMANT  (NT)  < FIMANT  or  (b)  is  a firm  track  with  time  since  last  update 
> TUPDT  and  XMANT  (NT)  < FMANT; 


• A track  is  limited  to  using,  three  regions  if  it  is  a firm  track  with  time  since 
last  update  > TUPDT  or  XMANT  (NT)  < FMANT; 


• A track  is  limited  to  using  two  regions  if  it  is  a tentative  track  with 
XMANT  (NT)  > FIMANT  or  is  a firm  track. 


A track  uses  the  largest  number  of  correlation  regions  for  which  it  qualifies.  The  con- 
ditions given  boil  down  to  the  conditions  that  the  track  has  been  moving  in  a straight 
line  for  some  time,  there  is  a good  velocity  measurement,  and  the  information  on  the 
track  is  current. 

One  last  function  of  the  correlation  routine  is  to  periodically  (every  TTR2CL  seconds) 
examine  the  average  range  velocity  of  each  track  to  see  if  it  should  be  converted  to  a 
clutter  point  (velocity  less  than  DIFTC). 


Track  Initiation 

The  method  of  track  initiation  is  the  weakest  area  of  this  particular  tracker.  This 
is  because  in  the  original  realtime  ADIT  system,  implemented  on  a Data  General  NOVA 
800  (a  minicomputer  with  an  SOO-ns  cycletime  and  32,000  16-bit  words),  there  was  in- 
sufficient time  and  space  to  implement  a more  sophisticated  initiation  algorithm.  However, 
because  an  investigation  of  initiation  techniques  is  one  of  the  first  items  to  be  attacked 
with  this  new  tool,  it  was  decided  to  adopt  the  ADIT  philosophy  with  minor  changes. 


The  initiation  technique  currently  used  is  as  follows:  Until  time  TICLUT,  only 
clutter  points  are  established  by  detections  left  in  the  input  file  after  the  correlation 
process  is  finished.  After  time  TICLUT  only  tentative  tracks  are  established.  At  the 
first  update  of  a tentative  track  more  than  TFIX  seconds  after  initiation,  a decision  is 
made,  based  on  track  tr^ectory  (velocity  and  height),  to  convert  the  tentative  track  to 
either  a firm  track  or  a clutter  point. 

In  an  environment  of  a low  density  of  detections  this  algorithm  is  highly  effective. 
Random  false  alarms  have  low  probability  of  correlating,  and  although  tentative  tracks 
are  initiated  for  each  false  alarm,  they  are  dropped  before  they  are  converted  to  firm 
tracks.  The  method  is  fast  and  simple  to  implement  and  there  is  no  storing  of  uncor- 
related detections,  no  generation  of  possible  tracks,  and  no  tasting  of  these  tracks  for 
reasonable  tr^ectories.  In  a scenario  involving  a raid  or  a region  of  high  false-alarm 
density  where  detections  may  have  multiple  correlations,  a more  sophisticated  technique 
involving  the  preceding  operations  will  be  necessary.  Such  a technique  will  be  discussed 
in  a later  section. 
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Output  Structurp  and  Sample  Output 


Curri'ntly  thero  is  no  intiTnal  evahtation  of  trark  quality.  All  evaluations  of  surveil- 
lance-radar detection  {wformance  and  tracking  prograiu  |)erformance  are  by  the  usi'r.  As 
an  aid  to  the  user  five  increasingly  detuiknl  levels  of  printout  are  available.  Level  0,  ilis- 
playiHl  in  Fig.  2a,  lists  each  opportunity  of  uixlating  a firm  track.  A zero  measured  range 
(third  column)  indicates  that,  for  whatever  reason,  the  track  was  not  u|Klateil.  U'vel  1, 
dis(>layed  in  Fig.  2b,  includes  the  information  of  level  0 plus  remarks  on  the  origins  of 
trucks  (i>romot«l  from  tentative  track,  transferrwl  from  clutter  file)  and  dispositions  of 
trucks  (droftiHx!  frt)m  file,  triuisferrtHl  to  clutter  file).  lA'vel  2,  displayed  in  Fig.  2c,  lists 
the  information  of  level  0 for  both  firm  ami  tentative  tracks  but  no  remarks.  Tentative 
tracks  are  idcntifiinl  by  the  addition  of  1000  to  the  track  number  (first  column).  Thus 
1005  is  tentative  truck  5;  if  the  truck  is  convertetl  to  a firm  track,  1005  iH'comes  5. 
lA'vel  3,  displayeti  in  Fig.  2d,  inclvides  the  level-2  updates  on  firm  and  tentative  tracks 
and,  in  analogy  to  level  1,  also  includes  remarks  on  the  origins  juid  disimsitions  of  the 
firm  and  tentative  tracks.  IawcI  4,  displayed  in  Fig.  2e,  eontains  the  information  of  level 
3 plus  information  on  clutter  points.  This  information  includes  the  creation  ami  deletitm 
of  clutter  points,  lists  all  detections  which  ixirrelate  with  the  clutter  point,  and  lists  the 
(>arameters  of  the  clutter  jmint  ami  the  detection  ustnl  to  ujxlate  it.  lA*vel  4,  when  ustxl 
on  data  derivwl  from  the  SURllKT  radar  (>rogram,  enables  the  usi'r  to  determine  what 
ha|ipens  to  every  iletection  presenttxl  to  tbt'  tracking  program.  This  is  (lossible  becau.se 
every  detection  from  the  SURDKT  program  is  ixkIihI  to  reveal  its  source:  valid  target, 
fixeil  clutter,  or  viuriable  clutter. 
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FirrURE  DEVELOPMENTS 
/Vnticlutter  and  ECM  Techniques 

In  an  ofH'rational  ADIT  s>’sU'm,  a desirable  feature  would  be  a methcxl  of  rontrolling 
clutter  or  ECM  false  alarms  in  the  detector  through  feedback  from  the  tracking  system. 

One  method  of  implementitig  this  would  be  to  maintain  a false-alarm  map  of  the  area 
covered  by  each  radar.  This  map  would  consist  of  the  numl>er  of  false  ahirms  in  each 
small  region  in  the  area.  The  numt)er  of  false  alarms  would  l>e  determinetl  by  tlie  tracking 
program,  jwssibly  simply  by  counting  the  number  of  detections  which  do  not  comdate 
with  any  track  or  fixed  clutter  fwint.  This  map  coulil  then  be  ustxi  to  raise  or  lower  the 
detection  thresholds  of  each  radar.  The  prt‘sent  program  is  not  rt‘al  time.  However,  be- 
cause part  of  the  detection  data  from  the  SlIRDET  program  is  the  signal-to-noise  ratio  of 
each  detection,  the  methixl  describetl  ixnild  be  evaluatixl  by  a postdettx'tion  filtering  of 
the  detection  file.  This  filtering  would  approximate  an  incn'astxl  thrt*sliold  by  deleting 
dettH'tiOns  with  low  signal-to-noist'  ratios. 

Initiation 

.•\s  discussixl  earlier,  a simple  initiation  scheme  can  dtyrade  the  track  file  in  a denst* 
detection  environment.  Thest'  iletections  ixiukl  be  f:Ust'  alarms  or  valid  detections  on  a 
multitarget  raid.  In  any  casi',  with  no  a-priori  knowUnlge  of  Uuget  motion,  huge  correla- 
tion regions  are  neivssar>’  to  allow  for  tiuget  motion.  This  allows  false  wnelations  to 
generate  tracks  on  random  falst'  ahums  or  to  confuse  targets  and  genende  spurious  vtdo- 
cities  in  a raid  environment. 

One  possible  candidate  system  would  ojx'rate  as  follows:  The  .system  would  l>e  a 
seixuate  subprogram,  divorml  from  the  track  correlation  and  iqxlating.  The  uiuxJirelatesl 
detections  remaining  in  the  files  after  the  wrrelation  pn>t'ess  would  be  saved  by  scan  and 
sector  in  a stqwate  file.  .\n  uiuxirndattnl  detection  in  the  current  scan  of  the  lu-ea  would 
trigger  a tnyectory-generation  event.  This  event  woukl  attempt  to  assix-iate  tletections 
from  the  last  n scans  with  a straight-line  trajectory.  (An  attempt  to  initiate  on  a maneu- 
wring  target  in  the  presetJtx*  of  fal.s«*  ahu'ms  and  misstxl  detections  would  rtxjuire  too  much 
time  real  time,  ni't  computation  time  -to  be  feasible.)  There  would  be  allowwl  at  most 
m missixl  tietections.  Because  of  the  tuHxl  to  identify  tracks  tpiickly.  n and  m would  be  a 
function  of  the  numlx'r  of  opportunities  to  detect  the  hyi>othetical  tnqectoiy  in  a givi'n 
amount  of  time.  One  implementation  of  this  process  would  Ix'  to  fit  a straight  line  to  each 
wllection  of  detections  (each  detection  from  the  appropriate  scan)  and  calculate  a rix>t- 
mean-s<(UBri'il  normaliztHi  error  (normalized  by  measurement  varianw.  averagtxl  over  the 
nun>lx?r  of  detectii>ns  us»’d)  and  conqwre  this  vidue  to  a tlm'shold.  The  wllection  with 
smallest  error  (less  than  the  threshold)  would  be  ilechutnl  a track.  The  detections  ustxl 
would  be  tleletixl  from  the  files.  As  the  numlx'r  of  detections  jx'r  scan  anil  the  numlx'r 
of  si'ans  increasi'.  (he  numlx'r  of  collections  inert'ast'  binomially  with  m and  n.  .A  sinqiler 
methixl  woukl  Ix'  to  fix  m aixl  n as  3 out  of  4 detections,  generate  straight  lines  from 
the  current  detection  to  each  detection  in  the  first  scan,  and  si'arch  the  iletections  in  the 
intervening  two  .si'ans  for  a closest  detection  which  is  sufficiently  dost'.  This  woukl  Ix' 
declanxl  a track,  anil  the  detections  would  be  removeil  from  the  files.  .\ny  detection 
more  than  4 scans  in  the  past  woukl  Ix'  iliscanltxl  as  a falsi'  alarm.  Due  to  the  large 
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measurement  variances  in  azimuth  and  elevation,  this  process  would  prolwbly  consider 
only  range  and  time  of  detection. 

An  example  of  the  gt^ometry  of  this  proi'ess  is  shown  in  Fig.  3.  where  detections  are 
located  in  a range-time  coonlinate  system  hy  a numlwr  indicating  in  which  scan  they 
occurred.  As  shown  in  the  diagram,  tr^ectory  A would  he  declared  a track  on  the  basis 
of  the  detections  of  scans  1,  2,  and  4,  but  an  extra  correlation  step  would  be  necessary 
to  eliminate  the  detection  of  scan  3,  which  evidently  is  a detection  of  this  target,  from 
the  files.  Triyectory  B would  be  declartxl  a track.  Trajectory  C does  not  have  an  inter- 
mediate detection  close  enough  to  declare  it  a track.  The  daslit^l  line  D would  not  even 
be  a candidate  trtyectory;  however,  if  it  is  detected  on  either  of  the  next  two  scans,  it 
would  then  be  declared  a track.  This  process  would  be  repeattxl  with  each  of  the  other 
detections  of  scan  4.  We  see  from  dashed  line  E that  possibly  an  inferior  track  has  been 
established  due  solely  to  the  order  in  which  the  detections  of  scan  4 were  processed.  To 
avoid  this  situation  would  greatly  complicate  the  process,  and  there  is  still  the  likelihood, 
as  in  case  D,  that  a track  would  still  be  established  on  the  next  two  si’ans.  The  problem 
of  the  possibly  erroneous  track  established  by  trjyectory  B ixiuld  be  solved  by  establish- 
ing tentative  tracks  from  these  tr^ectories  and  requiring  additional  confirmation  before 
promoting  these  tracks  to  firm  tracks. 


Fig.  3 - Example  of  initiation  geometry 


Correlation 

A problem  with  any  simple  wrrrelation  technique  is  that  of  wrong  associations  in  a 
dense  detection  environment.  If  the  correlation  technique  is  to  obtain  a track  and 
associate  the  track  with  the  nearest  detection,  a misseti  detection  can  cause  the  track  to 
be  u(Hlated  with  a detection  closer  to  another  track  and  in  fact  dorivetl  from  that  track. 
Conversly,  if  a detection  is  obtained  and  associattnl  with  the  nearest  track,  a false  alarm 
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(or  even  an  initial  detection  of  a new  track)  can  update  a track  which  has  been  detected  i 

and  whose  detection  is  closer  to  the  track  than  the  false  alarm.  A more  sophisticated 
correlation  scheme  is  desirable.  One  such  technique  has  been  described  in  the  section 
“Track  Correlation  and  Updating”  involving  multiple  nested  correlation  regions.  This  is 
a “discrete”  process  to  correlate  tracks  to  detections  and  then  detections  to  tracks.  This 
process  can  lead  to  enor  if,  for  example,  two  tracks  A and  B are  given  with  one  detection 
closer  to  B than  A and  track  A is  processed  first.  If  the  detection  lies  in  the  ith  correla- 
tion region  of  B and  in  the  (i  l)th  correlation  region  of  A,  then  the  detection  will  be 
flagged  as  being  closer  to  B and  track  A will  not  be  updated;  eventually  the  detection 
will  be  used  to  update  track  B.  If  the  detection  lies  in  the  ith  correlation  region  of  both 
A and  B,  then  the  first  track  considered.  A,  will  be  updated.  A correlation  scheme  used 
by  the  Applied  Physics  Laboratory  is  a “continuous”  process  in  which  a list  is  made  of 
tracks  due  to  be  updated  [9] . Associated  with  each  track  is  a list  of  detections  that 
correlate  with  the  track  and  the  corresponding  statistical  distance  between  each  detection 
and  track.  The  detection  with  smallest  statistical  distance  is  found  and  used  to  update 
the  corresponding  track.  That  track  and  detection  are  eliminated  from  the  lists,  and  the 
next-smallest  statistical  distance  is  found.  This  process  is  continued  until  either  the  tracks 
or  detections  are  exhausted.  One  modification  to  either  of  the  above  two  correlation 
schemes  is  simple  to  discriminate  against  tentative  tracks  in  favor  of  firm  tracks.  1 his 
modification  would  multiply  each  component  (range  or  azimuth)  of  the  correlation  error 
(the  difference  between  track  position  and  detection  position)  for  tentative  tracks  by  an 
“expansion”  factor  (such  as  1.26).  For  the  “continuous”  correlation,  detections  correlat- 
ing with  both  a firm  track  and  a tentative  track  are  more  likely  to  be  used  to  u(xlate  the 
firm  track.  For  the  "discrete”  correlator,  detections  whose  distance  from  a tentative 
track  is  more  than  80%  of  the  ith  correlation  size  would  not  correlate  until  the  (i  + l)th 
correlation  region. 

An  important  problem  as  yet  unaddressed  by  any  tracking  system  because  of  its 
complexity  is  illustrated  in  Fig.  4.  Given  tracks  A,  B,  and  C and  detections  1,  2,  and  3, 
the  list  technique  would  update  C with  2 and  B with  3,  leave  A unupdated,  and  interpret 
1 as  a false  alarm  or  initial  detection  of  a new  track.  The  question  is  whether  this  is  the 
most  reasonable  ass-  xsment  of  the  situation  or  whether  it  is  more  reasonable  to  assume 
that  A correlates  with  3,  B correlates  with  2,  and  C correlates  with  1 and  that  there  are 
no  false  alarms  or  missed  detections.  This  question  can  be  answered  using  maximum- 
likelihood  techniques  in  which  each  combination  of  tracks  and  detections  is  considered 
and  its  likelihood  calculated  based  on  the  probability  of  the  measurement  error  observed, 
probability  of  false  alarms  required,  probability  of  missing  detections,  probability  of  a 
new  track  being  detected  for  the  first  time,  and  even  an  estimate  of  the  probability  of  a 
target  maneuver  occurring  at  this  time.  This  method  of  solution  would  entail  much  more 
computing  time  than  is  normally  available  in  an  operational  system  but  may  be  of  interest 
as  an  evaluation  tool  in  determining  which  feasible  correlation  techniques  approach  this 
optimal  level. 

Filtering 

The  simplest  improvement  to  the  decreasing-bandwidth  filter  described  in  the  section 
“Tracking  Filter”  would  be  to  implement  the  adaptive-bandwidth  filter  used  in  an  earlier 
version  of  ADIT.  This  filter  was  not  discarded  for  any  shortcoming  but  in  the  interest 
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of  conserving  space  in  the  central  memory  of  the  minicomputer.  This  filter  is  described 
in  Ref.  4 and  in  greater  detail  in  Ref.  10.  The  principle  of  this  filter  is  to  set  the  filter 
bandwidth  by  the  magnitude  of  the  correlation  between  the  difference  between  measured 
and  pretlicted  (xjsition  on  the  current  ufxlate  and  the  difference  on  the  preceeding  ujxiate. 
If  the  track  carried  in  the  computer  closely  approximates  the  true  tr^ectory  of  the  target, 
then  the  measured  position  woukl  l>e  random  about  the  pretlicted  position,  the  correlation 
would  be  small,  and  the  bandwidth  should  be  correspondingly  small  and  smooth  the 
track  heavily.  If  the  track  is  not  the  true  triyectory,  due  to  a maneuver  for  example,  the 
measuretl  position  would  be  consistantly  displaced  relative  to  the  predictetl  position, 
causing  an  increased  correlation.  This  would  increase  the  bandwidth,  placing  more 
weight  on  the  measurements,  and  allow  the  filter  to  follow  the  target  through  the 
maneuver. 


Fig.  4 - Example  of  correlation  ambiguities 


Another  possibility  for  improving  the  filter  would  be  to  implement  a Kalman  filter 
with  either  a turn  detector  or  an  adaptive  bandwidth  as  described  above. 

Whatever  else  is  done  with  the  filter,  the  necessity  exists  for  tracking  in  a rectilinear 
(xy)  coordinate  system  at  short  ranges  (less  than  20  n.mi.). 

Software  Stabilization 

Naval  radars  generally  are  unstabilized;  the  antenna  shares  the  motion  of  the  plat- 
form. A 3D  radar  can  be  electronically  stabilized;  that  is,  if  the  location  of  the  target 
in  the  deckplane  coordinates  and  the  orientation  of  the  deckplane  are  known,  the  coor- 
dinates of  the  target  in  a stabilized  coordinate  system  can  be  calculated.  This  is  not 
possible  with  a 2D  radar,  because  the  elevation  of  the  target  is  not  known  and  the  mea- 
sured azimuth  of  the  target  is  a function  of  the  elevation.  However,  in  an  integrated 
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tracking  system  which  includes  2D  and  3D  radars,  the  elevation  information  can  be 
supplied  by  the  3D  radar  whenever  one  is  able  to  t'orrelate  the  unstabilizt'd  2D  detection 
with  a 3D  track.  One  will  then  l>e  able  to  stabilize  the  2D  information  iH'fore  using  it  to 
u|Kiate  the  track  or  even  before  correlation  is  attempted.  This  process  riHjuires  know- 
ledge of  the  platform  orientation.  This  information  is  available  in  data  generatetl  by  the 
SURDET  radar  simulation  model. 

An  even  more  interesting  problem  is  an  intt>gratetl  tracker  involving  only  unstabi- 
lized 2D  radars.  By  assuming  a straight  target  trajectory,  the  deviations  of  the  measureil 
azimuth  from  the  smoothed  azimuth  at  different  platform  orientations  would  yield  a 
gross  estimate  of  target  elevation. 

SUMMARY 

The  MERIT  system  as  currently  implemented  is  suitable  (in  coryunction  with  the 
SURDET  program)  for  evaluating  automatic  tracking  by  proposed  or  existing  radars.  The 
radars  can  be  considered  individually  or  as  members  of  a radar  suite.  Important  charac- 
teristics of  the  radars  are  accuracy  of  measurement,  probability  of  detection,  and  prob 
ability  of  false  alarm.  With  a-priori  knowletlge  of  the  target  sc-enario  and  environment, 
the  MERIT  system  can  examine  the  effect  of  these  radar  characteristics  on  a tracking 
system.  Effects  which  may  be  observed  are;  gross  track  inaccuracy  or  track  confusion 
due  to  poor  detection  accuracy,  loss  of  track  due  to  poor  probability  of  detection,  or 
saturation  of  the  system  or  generation  of  false  tracks  due  to  excessive  falsi'  alarms. 

The  MERIT  system  may  also  serve  as  an  investigative  tool  in  evaluating  the  cost  of 
various  improvements  and  the  effect  of  the  improvement  from  the  viewpoint  of  track 
quality.  The  improvements  currently  being  i-onsidered  are:  improved  initiation,  advanced 
correlation  techniques,  automatic  radar  control  in  adverse  environments,  optimal  realiz- 
able tracking  filters,  and  software  stabilization  of  2D  radars.  As  the  body  of  information 
increases  with  the  operation  of  MERIT  with  various  radar  systems  in  various  threat 
scanarios,  new  problems  may  be  identified  and  new  solutions  be  suggested  for  evaluation. 
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